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are reviewed and a method for the extrapolation of lattice data to the physical values of m u and rrid is outlined. 



1. INTRODUCTION 

The first crude estimates [Q for the magnitude 
of the three lightest quark masses appeared 25 
years ago (values in MeV): 
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Many papers dealing with the pattern of quark 
masses have been published since 

One of the reasons for being interested in an 
accurate determination of the mass values is that 
these are not understood at all - we need to know 
the numbers to test ideas that might lead to an 
understanding of the pattern, such as the rela- 
tions with the lepton masses that emerge from 
attempts at unifying the clcctrowcak and strong 
forces. Another reason is that the Standard 
Model must describe the low energy properties 
of the various particles to an amazing degree of 
accuracy. At low energies, the weak interactions 
are frozen and the neutrini decouple - the Stan- 
dard Model reduces to the gauge field theory of 
SU(3)xU(l). In the framework of this effective 
theory, the quark masses occur as parameters in 
the Lagrangian, together with the coupling con- 
stants e, g, the vacuum angle 6 and the masses of 
the charged leptons. This Lagrangian is supposed 
to describe the low energy structure of cold mat- 

J For an extensive list of refs., see Like most of these, 
the present brief report is based on the hypothesis that the 
quark condensate represents the leading order parameter 
of the spontaneously broken symmetry. The underlying 
physics issues are discussed in fe| and the refs. therein. 



ter to a very high degree of precision, provided 
the parameters occurring therein are accurately 
known. The size of the atoms, for instance, is 
determined by a B = Air / (e 2 m e ) and thus only in- 
volves parameters for which this is the case. 

In the following, I discuss the three light quark 
masses in terms of the magnitude of m s and of 
the ratios m u /md and m a jm^ which characterize 
their relative size. 



2. MAGNITUDE OF m s 

Apart from the lattice approach, the best de- 
terminations of the magnitude of m s rely on QCD 
sum rules jo . For a recent, comprehensive review, 
I refer to J7|. A detailed discussion of the method 
in application to the mass spectrum of the quarks 
was given nearly 20 years ago Q. The result for 
the MS running mass at scale n = 1 GeV obtained 
at that time was m s (l GeV) = 175 ± 55 MeV. 

The issue has been investigated in consider- 
able detail since then: Various versions of the 
sum rules for the pseudoscalar and scalar cor- 
relation functions were studied, as well as sum 
rules for two-point functions formed with bary- 
onic currents. More recently, an entirely inde- 
pendent sum rule determination of m s based on 
the strangeness content of the final state in r- 
decay became possible - in principle, this source 
of information is the cleanest one, because the 
hadronic input does not involve theoretical mod- 
els, but can be taken from experiment. The 
analysis of the r decay data described in H leads 
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toQ m s (2GeV) = 114 ± 23MeV. As discussed 
there and in ||, this result supersedes the earlier 
one in |jTo|| , which is based on the same experi- 
ment and on the same theoretical framework. 

The various sum rule results for m s are re- 
viewed in 0,0. In |Tl|j , the "world average" is 
quoted as m s (2GeV) = 118.9 ± 12.2 MeV, cor- 
roborating the above numbers with a significantly 
smaller error bar. It is notoriously difficult, how- 
ever, to account for the systematic errors of the 
various investigations |l^| - in my opinion, the 
number given does not fully cover these. 

The results based on the lattice approach, 
which were reported at this conference, cluster 
around m s (2 GeV) ~ 90 MeV, that is at the lower 
edge of the range obtained on the basis of sum 
rules. I think that the results are consistent with 
one another. In my opinion, this is important, 
because it shows that the magnitude of m a is un- 
derstood - we should not be content with numer- 
ical simulations of the theory, even if these will 
in the long run yield the most accurate results. 
A problem would arise, however, if the gradual 
reduction of the systematic errors to be attached 
to the lattice results should turn out to push the 
result downwards - that would be difficult to re- 
concile with the sum rules, in particular also in 
view of the lower bounds derived from the posi- 
tivity of the spectral functions Jl3| . 



3. PSEUDOSCALAR MASSES 

The best determinations of the relative size 
of m u , m d and m s rely on the fact that these 
masses happen to be small, so that the proper- 
ties of the theory may be analyzed by treating 
the quark mass term in the Hamiltonian of QCD 
as a perturbation. The Hamiltonian is split into 
two pieces: 

Hqpd — Hq + H i , 



2 Note the difference in scale: In recent times it has become 
customary to quote values at the scale /i = 2 GeV, because 
perturbation theory is under better control there. With 
the conversion factor used in H], the old value mentioned 
above corresponds to m a (2 GeV) = 125 ± 40 MeV. 



where Hq describes the three lightest quarks as 
massless and Hi is the corresponding mass term, 

Hi = J d 3 x{m u uu + m d 3d + m s ss} . 

Hq is invariant under the group SU(3) R xSU(3) L 
of independent flavour rotations of the right- and 
lefthanded quark fields. The symmetry is bro- 
ken spontaneously: The eigenstate of Hq with 
the lowest eigenvalue, |0), is invariant only un- 
der the subgroup SU(3) V C SU(3) n xSU(3) L . Ac- 
cordingly, the spectrum of Hq contains eight 
Goldstone bosons, tt ± , tt°, K ± , K°, K°, r]. The 
remaining levels form degenerate multiplets of 
SU(3) V of non-zero mass. 

The perturbation Hi splits the SU(3) mul- 
tiplets, in particular also the Goldstone boson 
octet. To first order, the squares of the meson 
masses are linear in m Ul m d ,m s and the symme- 
try fixes the coefficients up to a constant: 

Ml+ = (m u + m d )B + 0(m 2 ) , 
M%-+ = (m u + m s )BQ + 0{m 2 ) , 
Mf = (m d + m s )B Q + 0{m 2 ) . 

In the ratios M 2 + : ■ M^-o , the constant Bq 
drops out. Using the Dashen theorem to ac- 
count for the e.m. self energies (see section ||), 
these relations imply Q 
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4. SECOND ORDER MASS FORMULAE 

The leading order mass formulae are subject to 
corrections arising from contributions which are 
of second or higher order in the perturbation Hi . 
These can systematically be analyzed by means of 
the effective Lagrangian method JT^JT^] . In this 
approach, the quark and gluon fields of QCD are 
replaced by a set of pseudoscalar fields describing 
the degrees of freedom of the Goldstone bosons 
7r, K, rj. The effective Lagrangian only involves 
these fields and their derivatives, but contains an 
infinite string of vertices. For the calculation of 
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the pseudo-scalar masses to a given order in the 
perturbation Hi, however, only a finite subset 
contributes. The term Am, which describes the 
SU(3) corrections in the ratio Ml /Ml, 



-{l + A M + 0(m 2 )} , 



M K _ m^r 

M 2 m u + m d 

involves the two effective coupling constants L§, 
Lg, which occur in the derivative expansion of 
the effective Lagrangian at first non-leading order 



A M = 



8(M| - Ml) 
F 2 



(2L 8 - L 5 ) + xlogs . 



The term xl°S s stands for the logarithms char- 
acteristic of chiral perturbation theory. They 
arise because the spectrum of the unperturbed 
Hamiltonian Hq contains massless particles - the 
perturbation Hi generates infrared singularities. 
The coupling constant L$ also determines the 
SU(3) asymmetry in the decay constants, 



F K 4(M* - Ml) 

T = f 2 5 

f " TV 



Xlogs . 



The first order SU(3) correction in the mass ratio 
(Aff - M 2 K +)I{M 2 K - Ml) turns out to be the 
same as the one in M\fM 2 Q: 



M| 



M 2 K , 



m d 
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-{1 + A M + 0(m 2 )} 



Hence, the first order corrections drop out in the 
double ratio 



M 2 M 2 K0 



The observed values of the meson masses thus 
provide a tight constraint on one particular ratio 
of quark masses: 



Q 



-{l + 0(m 2 )} . 



(3) 



The constraint may be visualized by plotting the 
ratio rris/md versus m u /rrid [^7|. Dropping the 
higher order contributions, the resulting curve 
takes the form of an ellipse: 



1 



(4) 
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with Q as major semi-axis (the term rh 2 /m 2 s has 
been discarded, as it is numerically very small). 



5. VALUE OF Q 

The meson masses occurring in the double ra- 
tio (||) refer to pure QCD. The Dashen theorem 
states that in the chiral limit, the electromagnetic 
contributions to Ml + and to M 2 + are the same, 
while the self energies of K° and ir° vanish. Since 
the contribution to the mass difference between 
7T° and 7r + from md — m u is negligibly small, the 
masses in pure QCD are approximately given by 



(M 2 +) QCD = ( M 2 o) QCD =M 2 o 



(M- 
(M 2 K o 



,QCD 



M 2 K+ - 
M 2 K o , 



mL 
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where M^o , M^+ , M^o , M K + are the observed 
masses. Correcting for the electromagnetic self 
energies in this way, relation ^ yields Q = 24.2. 
For this value of the semi-axis, the ellipse passes 
through the point specified by Weinberg's mass 
ratios, eq. (|l|). The Dashen theorem is subject 
to corrections from higher order terms in the chi- 
ral expansion, which have been analyzed, but are 
not fully understood. Since the K° ~ K + mass 
difference is dominated by the contribution from 
md — m u , the uncertainties in the e.m. part do 
not very strongly affect the value of Q, but they 
do show up at the 5% level. For a detailed dis- 
cussion, I refer to [^8[ . 

The isospin violating decay r/ — > Sir allows one 
to measure the semi-axis in an entirely indepen- 
dent manner p9[ . The transition amplitude is 
much less sensitive to the uncertainties associ- 
ated with the electromagnetic interaction than 
the K° — K + mass difference: The e.m. contri- 
bution is suppressed by chiral symmetry and is 
negligibly small p0[| . The transition amplitude 
thus represents a sensitive probe of the symmetry 
breaking generated by md—m u . In fact, the decay 
rate is proportional to Q , with a factor that can 
be calculated up to corrections of second order in 
the quark masses. The result for the semi- axis 
is consistent with the information obtained from 
K° — K + but is more accurate [21 231: 



Q = 22.7±0.8 . 



(5) 
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6. KM AMBIGUITY 

Chiral perturbation theory thus fixes the 
isospin breaking parameter m u /m ( i in terms of 
the SU(3) breaking parameter m s /md, to within 
small uncertainties. The ratios themselves, that 
is the position on the ellipse, are a more subtle 
issue. Kaplan and Manohar fjl]] pointed out that 
the corrections to the lowest order result, eq. (Q), 
cannot be determined on purely phenomcnologi- 
cal grounds. They argued that these corrections 
might be large and that the w-quark might ac- 
tually be massless. A couple of years ago, this 
possibility was widely discussed in the literature 
|23|] , in view of the strong CP problem. 

The reason why phenomenology alone does not 
allow us to determine the two individual ratios 
beyond leading order is the following. The matrix 

to' = aim + a2(m + )~ 1 dct m 

transforms in the same manner as m. Symmetry 
does therefore not distinguish to' from m. For 
a real, diagonal mass matrix, the transformation 
law for m u , for instance, reads 

m' u = aim u +a 2 m d m s . (6) 

Since the effective theory exclusively exploits the 
symmetry properties of QCD, the above transfor- 
mation of the quark mass matrix does not change 
the form of the effective Lagrangian - the trans- 
formation may be absorbed in a suitable change of 
the effective coupling constants |17| ]. This implies, 
however, that the expressions for the masses of 
the pseudoscalars, for the scattering amplitudes 
or for the matrix elements of the vector and ax- 
ial currents, which follow from this Lagrangian, 
are invariant under the operation to — >m'. Con- 
versely, the experimental information on these ob- 
servables does not distinguish m u from m' u . In- 
deed, one readily checks that the transformation 
to — > to' maps the ellipse onto itself (up to terms 
of order {ra u — m^) 2 /to 2 , which were neglected). 
Since the position on the ellipse does not remain 
invariant, it cannot be extracted from these ob- 
servables within chiral perturbation theory. 

One is not dealing with a hidden symmetry of 
QCD here - this theory is not invariant under 
m — > m' . Some authors have made the claim that 



there is no reason for the quark masses entering 
the Lagrangian of the effective theory to be the 
same as the running masses of QCD. This claim is 
incorrect, as can explicitly be demonstrated with 
the following simple example. The Ward identity 
for the axial current implies that the vacuum-to- 
pion matrix element of the pseudoscalar density 
is given by^ 

(0\di~f 5 u\ir+) = V2F„Ml + /{m u + m d ) . 

The relation is exact, except for electroweak cor- 
rections (in the MS scheme, both the matrix ele- 
ment and the quark masses depend on the run- 
ning scale, but the relation holds at any scale). 
In the framework of the effective theory, the two 
sides can separately be calculated in terms of the 
effective coupling constants and the quark masses 
occurring in the effective Lagrangian. The rela- 
tion is obeyed, order by order in the chiral pertur- 
bation series. What appears on the right, how- 
ever, are the quark masses of the effective La- 
grangian. The effective theory thus reproduces 
the matrix element {Q\di^u\-K + ) if and only if 
the quark masses in the effective Lagrangian are 
identified with the running masses of QCD. 

7. THE RATIOS m u : m d AND m s : m d 

The KM-ambiguity is of phenomenological na- 
ture: Unfortunately, an experimental probe sen- 
sitive to the scalar or pseudoscalar currents does 
not exist - the electromagnetic and weak interac- 
tions happen to probe the low energy structure 
of QCD exclusively through vector and axial cur- 
rents. This means that theoretical input is needed 
to determine the size of the two individual mass 
ratios beyond leading order. 

One approach relies on the expansion in 1/N C . 
Since the KM-tranformation (^|) violates the 
Okubo-Iizuka-Zweig rule, the parameter a 2 is 
suppressed in the large N c limit. In order to 
cover this limit, the effective theory must be ex- 
tended, including the degrees of freedom of the r( 
among the dynamical variables, because this par- 
ticle becomes massless if the number of colours is 

3 The relation involves the matrix element of a pseu- 
doscalar operator — these are not KM-invariant, even at 
the level of the effective theory. 
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sent to infinity and the quark masses are turned 
off. In fact, within that framework, the KM- 
ambiguity disappears altogether: The transfor- 
mation m! — > to preserves the large N c counting 
rules only if «2 vanishes to all orders in \/N c [f24f . 
An evaluation of the mass ratio m s /rh based on 
the large N c analysis of the transitions rf — > 77 
and r] — ► 77 is described in pj|. 

Another method is based on the hypothesis 
that the lowest resonances dominate the low en- 
ergy behaviour of the various Green functions - 
once the poles and cuts due to the exchange of 
the Goldstone bosons are accounted for. More 
precisely, one first shows that the effective cou- 
pling constants can be represented in terms of 
convergent integrals over suitable spectral func- 
tions and then assumes that these sum rules are 
approximately saturated by resonances [^6| . 

Both of these methods lead to the conclusion 
that the corrections to the lowest order mass for- 
mulae are small: 

TJX TTi 

— = 0.553 ±0.043 , — = 18.9 ±0.8. (7) 

m d m d 

For a detailed discussion of this result, I refer to 
|l8[ . The values for the ratio (m u + md)/m s ob- 
tained by means of sum rules |ll[| as well as the 
lattice results presented at this conference are in 
remarkably good agreement with these numbers. 

8. EXTRAPOLATION IN m u AND m d 

The remainder of this report concerns the in- 
terface between the lattice approach and the ef- 
fective theory. As pointed out long ago (30|, chi- 
ral perturbation theory allows one to predict the 
dependence of the various correlation functions 
on the quark masses, as well as on the tempera- 
ture and on the volume of the box used. In the 
following, I only discuss the predictions for the 
dependence on the quark masses. 

By now, dynamical quarks with a mass of the 
order of the physical value of m s are within reach, 
but it is notoriously difficult to equip the two 
lightest quarks with their proper masses. Sup- 
pose that we keep m s fixed at the physical value 
and set m u = to^ — m, but vary the value of to 
in the range < to < i m s (at the upper end of 
that range, the pion mass is about 500 MeV). 



The expansions of M v and F n in powers of to 
are known to next-to-next-to leading order, from 
a two loop analysis |3l|-[33[ of the effective field 
theory based on SU(2) R xSU(2) L : 

M 2 

-JL = 1 -i x i 3+ rr. x Hl + x 2 k M + 0(x 3 ) , 
= 1 + xh - I x 2 ! 2 + x 2 k F + 0(x 3 ) , 

r 

l M = 77 ( 28 *i + 32 ^ -94 + 49) , (8) 
51 

h = —(144 + 164 + 64-64 + 23) , 

where M 2 = 2rhB is the leading termQ in the ex- 
pansion of M 2 , while F is the pion decay constant 
in the limit to = at fixed m s , and x denotes the 
dimensionless ratio 

In this representation, the chiral logarithms are 
hidden in the quantities £ n , which represent the 
running coupling constants at scale fi — M and 
depend logarithmically on M 

*n = ^§i- (10) 

The mass-independent terms k M and k F account 
for the remainder of 0(M 4 ), in particular for the 
contributions from the effective couplings of . 

The relations rt§) are dominated by the con- 
tributions of O(AP), which involve only two pa- 
rameters: the values of A3 and A4, respectively. 
In effect, these two parameters replace the coef- 
ficients a M , a F in the polynomial approximations 
Ml = M 2 (l + a M M 2 ), F n = F(l + a F M 2 ) that 
are sometimes used to perform the extrapolation 
of lattice data. In contrast to these approxima- 
tions, the formulae (||) do account for the relevant 
infrared singularities and are exact, up to and in- 
cluding 0(M 4 ). 

9. NUMERICAL DISCUSSION 

The estimates given in |Q confirm the expec- 
tation that the contributions from k M and k F are 

4 Note that B differs from the constant Bo occurring in 
section ts by terms of 0(m s ). 
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of order (M/M s ) 4 , where M s ~ 1 GeV is the 
mass scale characteristic of the scalar or pseu- 
doscalar non-Goldstone states contributing to the 
relevant sum rules (in the SU(2) framework we are 
using here, the KK continuum also contributes 
to the effective coupling constants, but in view of 
4M|- ~ M|, the relevant scale is even somewhat 
larger) . Unless the quark masses are taken much 
larger than in nature, these terms are very small 
and can just as well be dropped. The relations 
(g) then specify M„ and F„ as functions of rh, in 
terms of the 6 constants F, B, Ax, . . . , A4. 

Except for B and A3, all of these can be de- 
termined quite accurately from experiment: The 
scales Ax and A 2 can be evaluated on the ba- 
sis of the Roy equations for tttt scattering, while 
A4 is related to the scalar charge radius, which 
can be extracted from a dispersive analysis of 
the scalar form factor j3^|5(|. For the purpose 
of illustration, I use the lowest order relation 
B ~ {Mjc — \M*)/m B to replace the variable 
rh by the ratio rh/m s . For A 3 , I invoke the crude 
estimate lnA|/M| = 2. 9 ±2. 4 given in [§6), which 
amounts to 0.2 GeV < A 3 < 2 GeV. 

Consider hrst the ratio F n /F, for which the 
poorly known scale A3 only enters at next-to- 
next-to leading order. The upper one of the two 
shaded regions in fig. 1 shows the behaviour of 
this ratio as a function of rh, according to for- 
mula (||). The change in F^ occurring if rh is 
increased from the physical value to \ m s is of 
the expected size, comparable to the difference 
between Fk and F„. The curvature makes it ev- 
ident that a linear extrapolation in rh is mean- 
ingless. The essential parameter here is the scale 
A4 that determines the magnitude of the term of 
order M 2 . The corrections of order M 4 are small 
- the scale relevant for these is Af — 0.5 GeV. 

In the case of the ratio M 2 /M 2 , on the other 
hand, the dominating contribution is determined 
by the scale A3. The fact that the informa- 
tion about that scale is very meagre shows up 
through very large uncertainties. In particular, 
with A 3 ~ 0.5 GeV, the ratio M 2 /M 2 would re- 
main very close to 1, on the entire interval shown. 
The corrections of 0(M 4 ) are small also in this 
case (the relevant scale is Am — 0.6 GeV). 

The above discussion shows that brute force is 
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Figure 1. Dependence of the ratios F^/F and 
M 2 /M 2 on rh — \{m u + mj. The strange quark 
mass is held fixed at the physical value. The ver- 
tical line corresponds to the physical value of rh. 

not the only way to reach the very small values 
of m u and rrid observed in nature. It suffices to 
equip the strange quark with the physical value 
of m s and to measure the dependence of the pion 
mass on m u ,md in the region where M n takes a 
value like 400 or 500 MeV. Since the dependence 
on the quark masses is known rather accurately 
in terms of the two constants B and A3, a fit 
to the data based on eq. (||) should provide an 
extrapolation to the physical quark masses that 
is under good control. Moreover, the resulting 
value for A3 would be of considerable interest, 
because that scale also shows up in other con- 
texts, in the tttt scattering lengths, for example. 
A measurement of the mass dependence of F w in 
the same region would be useful too, because it 
would provide a check on the dispersive analysis 
of the scalar radius that underlies the determia- 
tion of A4 (in view of the strong unitarity cut in 
the scalar form factor, a direct evaluation of the 
scalar radius on the lattice is likely more difficult). 

10. 7T7T SCATTERING LENGTHS 

As a further example, I consider the two tttt 
5-wave scattering lengths, for which the depen- 
dence on the quark masses was discussed already 
in J37J , including a comparison with the quenched 
lattice data available at the time. The two loop 
representation |Q explicitly specifies the scat- 
tering lengths in terms of the effective coupling 
constants, up to and including contributions of 
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Figure 2. 7T7t S'-wave scattering lengths as a func- 
tion of the pion mass. 

0(Ajf 6 ). Fig. 2 shows the corresponding correc- 
tions to Weinberg's low energy theorem^ fl38|] 
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1 = 2 



Ml 
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Note that the Weinberg formulae are written in 
terms of the values of i 7 ^ and that corre- 
spond to the quark mass of interest - instead of 
the corresponding lowest order terms F and M. 
Although the range shown here is considerably 
smaller than in fig. 1, the corrections are much 
larger: The expansion of the scattering lengths 
in powers of M n converges only very slowly. The 
effect arises from the chiral logarithms associated 
with the unitarity cut, which in the case of the 
1 = channel happen to pick up large coefficients. 
In fact, the chiral perturbation theory formulae 
underlying the figure are meaningful only in the 
range where the corrections are small (the shaded 
regions exclusively account for the uncertainties 
in the values of the coupling constants). 

I conclude that, in the case of the 1 = scat- 
tering length, the extrapolation in m requires sig- 
nificantly smaller quark masses than the one for 
Mir or F^. In the 1 = 2 channel, the effects are 
much smaller, because this channel is exotic: The 
final state interaction is weak and repulsive. 

The method proposed in |36| replaces the ex- 
pansion of the scattering amplitude at threshold 

5 The standard definition of the scattering length corre- 
sponds to ao/M n . It is not suitable, because it differs 
from the invariant scattering amplitude at threshold by a 
kinematic factor that diverges in the chiral limit. 



by one in the unphysical region, where the con- 
vergence properties are similar to those for F n 
and M%. Indeed, that method yields a remark- 
ably precise prediction for the scattering lengths: 
a% = 0.220 ± 0.005, = -0.0444 ± 0.0010. The 
lattice result, ag = -0.0374±0.0049 (39), is on the 
low side, but not inconsistent with the prediction. 

11. CONCLUSIONS 

The results for m s obtained from recent lattice 
simulations with dynamical quarks are consistent 
with those based on sum rules, but the central 
values are on the low side. 

The lattice results for the ratios m u /md and 
Tn s /md confirm the values obtained in the frame- 
work of chiral perturbation theory. 

The dependence of the various observables on 
the quark masses contains infrared singularities 
which can be worked out by means of chiral 
perturbation theory. This information might be 
useful for the extrapolation of lattice results to 
the small quark mass values of physical interest. 
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of Science, G. Colangelo, J. Gasser and S. Nar- 
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